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Abstract 

The multiplicity distributions for individual fragment Z values in nuclear 
multifragmentation are binomial. The extracted maximum value of the mul- 
tiplicity, mzi is found to depend on Z according to mz = Zq/Z, where Zq is 
the source size. This is shown to be a strong indication of statistical coverage 
of fragmentation phase space. The inferred source sizes coincide with those 
extracted from the analysis of fixed multiplicity charge distributions. 
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Is nuclear multifragmentation |I]j2| a statistical or a dynamical process? Let us consider 
a fragmentation space defining the number and masses/charges of the produced fragments. 
The thermal population of each cell can be expressed in terms of suitable Boltzmann factors. 
The thermal nature of this population can be tested in a variety of ways. A very visual way is 
the Arrhenius plot where the log of the population is plotted vs. the reciprocal temperature 



The question however, remains whether, apart from the Boltzmann factor, all the cells 
of the fragmentation space are uniformly explored. A test of such a uniform filling implies 
the knowledge of the size (mass/charge) of the source. A new empirical feature observed in 
many reactions has led us to a way to verify the uniform filling of fragmentation space and 
to determine simultaneously the source size. 

It has been shown [|3| that the intermediate mass fragment (IMF) multiplicity distribution 
P n at any given transverse energy E t is empirically given by a binomial distribution 



This implies that fragments are emitted nearly independently of each other, so that the 
probability P n of observing n fragments can be written by combining a single one-fragment 



represents the maximum possible number of fragments, which is immediately related to the 
source size. 

The simplest statistical equilibrium model of multifragmentation has exactly the struc- 
ture of Eq. ([]]). Let us assume that the source is made up of m fragments. The "outside" 
fragments have energy €2, and those "inside" have energy e\. A generic partition of n frag- 
ments outside and m — n fragments inside has the probability: 



l/T. 
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emission probability p according to Eq. ([!]). The parameter m (the total number of throws) 
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which leads to Eq. ([!]) when 
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Thus, a simple way to obtain the size of the source is to multiply m by the fragment size. 

When the definition of IMF covers a range of atomic numbers (IMF: Zth < Z < 20, with 
Zth equal to 3), one should multiply m by a "suitably" averaged Z. In fact, a dependence 
of m on the lower threshold Z t h has been found such that m{Z t h) x Z t h ~ constant ||. The 
natural next step is to restrict the fragment definition to a single atomic number Z. 

A straightforward generalization of Eq. (p]) to the production of fragments with charges 
1, 2, ...Zq is given by the multinomial distribution 

p = ^r — \vrv7-v n z z : (4) 

n l \n 2 \...nz \ 

with 

Z = ^2 Zn z . (5) 

z 

Here there is no single quantity m as in Eq. ([]]), since the constraint is now on the total 
charge rather than on the total number of fragments. However, this leads immediately to 
a simple scaling that must be obeyed if the fragmentation phase space is to be completely 
explored. 

Let us first consider, as a familiar example, the fragmentation space spanned by the Euler 
number partitions. This case is particularly simple because all partitions are unbiased, or 
equally probable. This would correspond, for the different model defined in Eqs. (|j) and 
(§), to setting pj=constant. 

Let us indicate with W(Zq) the statistical weight (number of partitions) associated with 
the integer or "charge" Zq. If we now select the partitions containing at least nz integers of 
size Z, their number is given by the number of partitions of the residue 

P nz = W(Z - n z Z) (6) 

by definition. Since the "charge" constraint is applied "minimally," what counts is only the 
product Zn Zl rather than the individual Z values. Consequently, the weight of the residue 
does not change if we substitute n z Z with nz>Z' provided that n z Z = nz>Z'. In other 
words, P nz>z = Pn z ,,z> \tn z Z = n z >Z' '. 
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This gives immediately the scaling laws 



n z /n z > = Z' 



/z 



(7) 



or for the extreme value of nz = mz. 



m z = Zq/Z. 



(8) 



This result follows directly from the uniform exploration of the fragmentation phase space. 
It amounts to a complete decoupling of all the fragmentation paths. Once the system attains 
the charge loss of q = nZ, it is indifferent to how this was attained. Any path q = J2 niZi 
is equivalent and substitutable. Similarly, for the remaining charge Z — q. 

This argument works when all partitions are unbiased, as for instance in the Euler number 
partition mentioned above. However, in the binomial/multinomial distributions of Eqs. ([I]) 
and (H), each partition is weighted by the probabilities pi (reflecting Q- value effects) which 
may or may not be Boltzmann factors. 

Yet, the binomial/multinomial analysis elegantly and automatically separates out the Q- 
value dependent probability p from the parameter mz which now contains direct information 
on the accessible fragmentation space and must follow the scaling given by Eq. (|^). Thus the 
1/Z scaling is the counterpart of the Arrhenius plot which verifies the thermal (Canonical) 
population of each fragmentation configuration. 

When the restriction to individual Z values is made experimentally jp, the multiplicity 
distributions are found to be nearly Poissonian, namely mp « m. This introduces inter- 
esting simplifications in the analysis and interpretation of the data, but at the cost of a 
loss of scale. In the Poisson limit the average multiplicity (n) = mp is the only accessible 
parameter, and the decomposition into m and p becomes impossible. The recovery of scale 
for an individual Z is highly desirable in view of the possibility that the number of throws 
mz (for a single species Z) might obey the simple scaling of Eq. (^) . 

Thus, we have attempted binomial fits of the multiplicity distributions for individual 
Z values in an effort to extract mz- Fortunately, a number of reactions ( 36 Ar+ 197 Au at 
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35 to 110 AMeV §, 129 Xe+ 27 Al, 51 V, nat Cu, 89 Y, at 50 AMeV §, and 129 Xe+ 197 Au at 50 
to 60 AMeV J7|) have been studied with good Z resolution and high statistics. We first 
consider the asymmetric, intermediate-energy reactions in reverse kinematics exemplified by 
129 Xe+ nat Cu at 50 AMeV, for which we can expect a single dominant fragment source. 

Examples of both binomial and Poisson fits to the carbon yield from this reaction are 
shown in panel a) of Fig. |l[ An improvement of the fit by using the binomial expression 
is observed for large fold numbers. A similar improvement is observed for each Z in all 
reactions listed in this letter. 

The E t dependence of the parameters mz from the binomial fits to the multiplicity 
distributions associated with each fragment atomic number leads to several observations. For 
each Z value, mz increases to a near constant value with increasing E t . We approximate this 
behavior with the form mz = m^tanh/i^. The parameter m° z represents the saturation 
value of mz for large E t and / controls the rise of mz with increasing E t . The solid lines 
in panel b) of Fig. [I] are the empirical fits to mz values extracted for lithium and oxygen 
emission from the reaction 129 Xe+ nat Cu at 50 AMeV. The other discontinuous lines are fits 
to data not shown (Z=4-7). 

Furthermore, at all E t values there is an overall decrease of mz with increasing fragment 
Z value in agreement with the expected scaling Zmz = Zq. This remarkable dependence is 
exemplified in panel c) of Fig. [I] and in Fig. |2|. By applying the expected scaling (Zmz), all 
of the fits to the 129 Xe+ nat Cu data collapse together, resulting in the approximate source 
"size" as a function of E t . A weighted average {(Z )) of the data over different exit channels, 
constructed according to 

(Z }(E t )=Y,Zm z (E t )a z , (9) 

z 

is shown by the symbols in panel c) of Fig. [I]. The weight az is the standard weight 
(proportional to the inverse square of the individual errors). A similar behavior is observed 
in two additional asymmetric reactions 129 Xe+ 51 V and 89 Y (see Fig. |T]d). 

The E t dependence of the source size is tantalizing. The source size increases quickly 
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to a saturation value. The fact that E t is related to impact parameter as well as to the 
total excitation energy may explain the observed features. In the highly asymmetric reverse 
kinematic reactions one quickly achieves sufficient overlap to produce a dominant Xe-like 
source as one moves from peripheral to central collisions. 

As a special case of the 1/Z scaling, the "saturation" mz values from central collisions 
of the reverse kinematics reactions (the top 5% of the E t scale, shown by the hatched 
regions in panel b) of Fig. [I]), are shown in Fig. ^|. The open symbols represent the scaled 
quantity Zmz = Zq. The solid lines are weighted averages for the different reactions. 
The same data in the form tuz/Zq vs. Z are shown by the solid symbols, where the 1/Z 
dependence is manifested by the good agreement of the data with the values of 1/Z (solid 
line). This striking 1/Z dependence of the parameter mz is observed for all asymmetric 
systems measured. 

These overall results for asymmetric reactions suggest the dominance of a single source, 
strongly support the hypothesis of uniform (statistical) exploration of the fragmentation phase 
space, and lead to the interpretation of Z = Zmz as the source "size." The 1/Z scaling is 
general and should be observed in models used to describe multifragmentation. 

In the less asymmetric reactions 129 Xe+ 197 Au at 50 and 60 AMeV for which at least two 
sources are plausible, we shall refer directly to Z = Zmz as the source size, although we 
shall see that now it depends on the fragment Z value as well as on E t . The weak decrease 
of the source size with increasing fragment size Z, already observable in 129 Xe+ 89 Y (Fig. 0), 
becomes more visible in the case of the 197 Au target. At low fragment Z values, the source 
size Zq is ~ 70 and it decreases monotonically with increasing fragment size Z to to a source 
size Zq of approximately 40-50. This fragment size dependence seems to suggest that for 
the 89 Y target and, most of all, for the 197 Au target there may be a distribution of sizes, the 
higher Z fragments being emitted preferentially by the smaller source(s). 

The reactions 36 Ar+ 197 Au at 35, 50, 80, 110 AMeV give a picture intermediate between 
the 129 Xe+ 197 Au and the 129 Xe induced reverse kinematics reactions. They also give infor- 
mation of the source size dependence on bombarding energy. The source size at low fragment 
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Z increases from Zq « 30 to Zq m 60 as the bombarding energy increases, AMeV, while at 
higher fragment Z the source size increases from Z w 20 to Z w 40. 

In previous work ||, it was empirically shown that the observed charge distributions 
resulting from nuclear multifragmentation obey the following invariant form: 



where n is the total intermediate mass fragment multiplicity of the event; E t is the total 
transverse energy; and B(Z) is the "barrier" distribution. 

From thermodynamic considerations and percolation simulations, it was shown that c in 
Eq. (|T0| ) vanishes when the gas of IMFs is in equilibrium with a liquid (residue source, or 
percolating cluster), and assumes a value oc 1/Z (Z being the source size) when the source 
is wholly vaporized |||8|||. 

Experimentally, the parameter c undergoes an evolution with (transverse) energy from 
approximately zero to a positive constant || . Thus the source size evolves from near infinity 
(an "infinite" reservoir of fragments) to the actual size of the source. With the exception of 
the reactions 129 Xe+ 27 Al, 51 V, Cu, 89 Y, this limit is attained. Thus, it is possible to plot the 
value of Z determined from mz against that obtained from the c parameter (both for the 
top 5% most central collisions in E t ). 

Such a plot is shown in Fig. |j. The result is striking. Not only are the two quantities 
well correlated, but they also agree quite well in absolute value. This good agreement gives 
confidence that we have gained direct access to the source size. This source (sources) is 
specifically the entity that generates the fragments through "chemical equilibrium" . It does 
not contain the pre-equilibrium part which is often incorporated in other source reconstruc- 
tion methods. 

In conclusion, we have shown that: a) the binomial (nearly Poissonian) multiplicity 
distributions for individual fragment atomic numbers permit the extraction of the parameter 
mz, the number of throws; b) mz for reactions where a single source is clearly dominant 
has the form mz = Zq/Z; c) the 1/Z dependence is dramatic proof that the fragmentation 




(10) 
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phase space is statistically explored; d) source size(s) can be extracted and should reflect 
the region(s) where chemical (as opposed to physical) equilibrium is achieved; e) these 
source sizes agree with the sizes obtained from the analysis of multiplicity selected charge 
distributions in the E t range where a single gas phase, or thermodynamic bivariance, prevails. 
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FIGURES 
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FIG. 1. Panel a): the n-fold probability distribution (symbols) and Poisson (dashed) and bi- 
nomial (solid) fits are plotted as a function of transverse energy for carbon fragments emitted 
from the reaction 129 Xe+ nat Cu at 50 AMeV. Panel b): the extracted binomial parameter mz (the 
number of "throws") is plotted as a function of transverse energy for lithium (solid circles) and 
oxygen (open circles) emission. The solid lines are hyperbolic tangent fits to the indicated data. 
The other lines are fits to data not shown. The two hatched regions represent weighted averages 
of the top 5% most central collisions (based on the integrated Et spectrum) of the mz values for 
lithium and oxygen. Panel c): The fits from panel b) are scaled by the atomic number Z of the 
emitted particle. The square symbols represent an "average" source size calculated with Eq. @. 
Panel d): The symbols represent an "average" source size calculated with Eq. (Q) for 129 Xe+ 51 V 
and 129 Xe+ 89 Y. The lines are hyperbolic tangent fits to the data. 
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FIG. 2. Top panel: The mz values scaled by Z as a function of Z for the indicated reactions 
(open symbols) at 50 ^4MeV extracted for the 5% most central collisions. Bottom panel: The 
values of mz divided by the value of the extracted source charge Zq for the indicated reactions 
(solid symbols). The solid line is 1/Z. 
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FIG. 3. The scaled values Z x are plotted for the indicated reactions of the 5% most central 
collisions. The solid lines represent weighted averages from Z = 3 to 6. 
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FIG. 4. The source size Zo extracted from the binomial analysis as a function of the source size 
(1/c) determined as per ref. || from central collision of the indicated reactions. 
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